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16. Abstract 


A technique consisting of a marriage between stress freezing photo-elasticity 
and a numerical method was used to obtain stress intensity factors for natural 


cracks emanating from the corner at which a hole intersects a plate surface. 
Geometries studied were: crack depth to thickness ratios of approximately 0.2, 

0.5 and 0.75; crack depth to crack length ratios of approximately 1.0 to 2.0; and 
crack length to hole radius ratios of about 0.5 to 2.0. All final crack geometries 
were grown under monotonic loading and growth was not self similar, with most of 
the growth occuring through the thickness under remote extension. Stress intensity 
factors were determined at the intersection of the flaw border (i) with the 


plate surface (K_) and (ii) with the edge of the hole (K H ). Results showed that 
for the relatively shallow flaws K H = 1.51C, for the moderately deep flaws K„ = K $ 
and for the deep flaws K H * 0.5 1C, revealing a severe sensitivity of K to flaw 
geometry. Results were comparedVith the Bowie Theory, and approximate criteria 
developed by Hall and Finger and Hsu and Liu. These comparisons showed that these 
theories significantly over-estimated the SIF for moderately deep flaws a/t ® 0.5 
at both the plate surface and the hole, but for shallow flaws, and Hall-Finger theory 
under-estimated the SIF at the hole and the Bowie theory under-estimated the SIF at 


the surface for deep flaws. 
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- Stress components in plane normal to crack border (psi) 

- Normal stress in direction of crack extension near crack 
tip (psi) 

- Remote stress (psi) 

- Mode I stress intensity factor (psi - [in] ' ) 

- Polar coordinates (in, radians) 

- Flaw depth (in) 

- Flaw length (in) 

- Hole radius (in) 

- Plate thickness (in) 

- Fringe order 

- Material fringe value (lbs/in/order) (f = 1.46) 

- Mode I SIF at plate surface (psi - [in]^ 2 ) 

- Mode I SIF at hole surface (psi - [in] 1 / 2 ) 

- Mode I SIF Experimental (psi - [in] 1 / 2 ) 

- Mode I SIF Apparent (psi - [in] 1 / 2 ) 
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INTRODUCTION 


A common cracked body problem in the aerospace industry consists of a 
corner crack emanating from the intersection edge between a plate surface 
and a hole. No analytical solution exists for this problem even as a 
near field solution only. Yet the designer is forced to design against this 
type of crack using very approximate methods [1]. Apparently, the first 
study undertaken of the problem was an experimental study by Hall and 
Finger [2]. They inserted artificial flaws by EDM with depth/length 
less than unity and used fatigue loads to initiate cracks after which 
residual static strength tests were run. In evaluating their results, 
they assumed that a state of plane strain existed at the point of intersection 
of the flaw border and the boundary, but in their fracture criterion they 
had only one value for the stress intensity factor (SIF) and did not 
account for the variation of the SIF and constraint along the flaw border. 

As an alternate approach, they suggested modelling the corner flaw with an 
"equivalent" Bowie type [3] through crack. The latter approach was refined 
by Liu [4] but still required the selection of an arbitrary equivalent 
crack length for the Bowie model. This approach has recently been expanded 
by Hsu and Liu [5], and while this latter analysis contains some two 
dimensional approximations, and takes no account of back surface effects, 
it constitutes the most recent effort at quantifying the problem. 

Because of the technological importance of the above problem, and the 
approximate nature of existing analytical approaches necessitated by the 
intractability of the mathematical model of the complete problem, one is 
led to consider experimental techniques as an alternate approach. Stress 
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freezing photoelasticity is a well known technique for evaluating three 
dimensional stress distributions in the vicinity of stress raisers. A 
method for extracting the stress intensity factor from photoelastic data 
for two dimensional problems was proposed by Irwin [6] in 1953, and the 
method has been modified and refined substantially since that time. Recent 
studies by Kobayashi and his associates [7]— [10] have applied the method 
to dynamic photoelasticity, and the senior author and his associates 
[1 1 ]-[21 ] have refined the method for three dimensional problems. 

It was the purpose of the present study to determine stress intensity 
factors at the end points of flaws emanating from the corner formed by 
the intersection of a plate with a hole using stress freezing photoelas- 
ticity and to compare the results with the studies noted above. The authors 
used a numerical technique known as the Taylor Series Correction Method 
(TSCM) [20] in order to extract the SIF values from the photoelastic data. 
Before describing the experiments, a brief review of the analytical back- 
ground appears to be desirable. 


ANALYTICAL CONSIDERATIONS 

Consider the Irwin two parameter near field equations for Mode I 
loading: 

c = — — 1/2 cos -j (1 - sin|- sin ^ ) - a 
nn ( 2 »r )' /2 2 2 2 0 

K 

o = — — cos *■ (1 + sin sin ^ ) (1) 

“ <2*r) ,/2 

T = — — , /? S1 ’ n f COS | C0S^| 

nz <2ur) 1/2 2 2 2 
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where the notation, adapted to the problem at hand, is pictured in Figure 
1. Here o 0 is the part of the regular stress field which is independent 
of r. Although Equations (1) were originally proposed for two dimensional 
problems, the singular parts have been shown to be valid for stresses in 
planes perpendicular to el 1 i pti cal ly shaped crack borders [22] and are 
generally applied to arbitrary shaped plane crack borders as well. By 
substituting Equations (1) into Equation (2): 


Sax' * ««m, ’ «zi>* * O'" 


nz 


and evaluating t v along 9 = tt/ 2, one obtains: 

iTlaA 


Ky K T o a /9 

— r I i I o i /o j_ o l * / 

T max - ^ 8irr 4(-rrr) ' 4 


( 2 ) 


(3) 


which may be combined with the stress optic law: 

_ Nf 
T max ” 2t 

in order to obtain an expression of the form: 


(4) 


Kj = f(cr 


O- N i- r 1> 


(5) 


from which Kj may be evaluated from experimental data N., r^ along 
9 = tt/2 (See references [16] and [17]) The stresses are evaluated along 
9 = tt/ 2 since the fringes spread out in approximately this direction and can 
be most accurately discriminated along this line. An alternate form of 
Equation (3) may be written as: 


max 


- A ‘ + 
' 7/2 + 


B'r 1 ^ 


( 6 ) 
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where 


A' = — and B-J = B' (A’, B') 

/§tT 0 

When one accumulates photoelastic data along e = tt/2, if the foregoing 
theory is to hold, the stresses in the data zone must be dominated by the 
singular stresses as given in Equations (1). Experience with three 
dimensional problems has revealed that the singular zone is often severely 
constricted and a part of the data may lie outside that zone. In order to 
adequately account for this effect, the authors have used a Taylor Series 
Expansion of the regular part of the maximum in plane shearing stress in the 
form: 


T max 




(7) 


where A = KjMk as before. 

If we ignore the dependence between the coefficients A', B^, B-j 
in Equation (6) we see that Equation (7) reduces to Equation (6) for a 
two degree of freedom system of equations to within truncation error. 
Moreover, Equation (7) corresponds to an application of the Williams Stress 
Function along 0 = tt/2 for two dimensional problems. 

In order to apply Equation (7), one determines the coefficients 
A, B n from a least squares analysis of the experimental data using a 
truncated form of Equation (7). Normally the lowest order curve which best 
fits the experimental data is used. Details are described in references 
[17] and [21]. 
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THE EXPERIMENTS 


Models 

A series of stress freezing photoelastic experiments were designed for 
the purpose of obtaining estimates of the SIF at points S and H along the 
flaw border. (Figure 1) The specimens were made from PLM-4B, a stress 
freezing material manufactured by Photolastic, Inc., Malvern, Pa. and Hysol 4290 
made by Hysol Corp. , Olean, N. Y. using the following procedure: 

i) Mill 30 mils from both surfaces of the plates, 

ii) Drill and ream holes. 

iii) Mount a razor blade in a special jig and tap in a quarter circular 
crack. 

iv) Mount test specimen in a stress freezing oven in a dead loading 
rig (Figure 2 ) and heat to critical temperature, 

v) Load with enough load to slowly extend flaw to its desired 
depth. Remove the load. 

i Test geometries studied are found in the upper part of Table I. 

It was originally intended to attempt to duplicate the geometries of the 
Hall and Finger tests. However, since natural cracks were used, only one 
crack dimension could be controlled, and the other dimension grew to its 
"natural" companion value. This led to flaws for which a/c > 1 .0 and self 
similar flaw extension did not occur. Instead, most of the growth occurred 
in the depth direction. In one test, (Test 9) a crack with a/c < 1.0 
was produced by flexing the plate to enlarge the c dimension. 

Test Procedure 

After crack growth was completed, live loads below the threshold value 
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to cause crack growth were applied above the critical temperature, 
and specimens were cooled under load, thus freezing in the fringes and 
deformations. Slices were then taken i) parallel to the plate surface for 
use in determining the SIF at point S and ii) tangent to the hole for 
determining the SIF at point H. The latter slice was sanded on the hole 
surface to constant thickness. Slice locations are shown in Figure 3. 

These slices were placed in a tank of liquid of the same index of refraction 
as the model material and fringe patterns were obtained using partial mirror 
fringe multiplication and also by the Tardy method. A typical fringe pattern 
for the hole H is shown in Figure 4. Alices were 0.02 to 0.07 in. thick. 

Results 

A typical set of data for points S and H are shown in Figure 5 
together with the TSCM curves generated from the data. These curves may 
be regarded as typical of all of the tests analyzed. That is, for all 
tests, the surface slices revealed linear data which indicates a two 

j 

degree of freedom system of equations while the slices along the hole 
revealed non-linear data. For these cases, however, the Tardy method 
revealed that this data was also linear over a portion of the data zone 
nearer to the crack tip. Thus, instead of using a higher order curve 
to fit all of the data, the authors elected to use a two degree of freedom 
system on that portion of the data in the linear range as shown for 
the H curve in Figure 5. The authors interpret the linear part of the 
curve to represent the singular zone and believe that extrapolation of 
the linear curve results in much less error than would be the case for a 
higher order curve. In order to verify this technique, a test was run 
where the corner crack was allowed to grow through the thickness of the 
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plate so as to form a Bowie type crack on one side of the hole. The hole 

radius was 0.19 inches and the crack length was 1.15 inches. The resulting 

K was 2 % higher than the result predicted by the Bowie analysis. The 
Ex 

authors estimate, however, that experimental results generally can vary 
by 5%. Even after Tardy analysis, the slopes of the curves of normalized 
K. vs. (-)^ 2 or were, in some cases, sufficiently large to cause 

the authors some concern as to the accuracy of the K^ £x values. In 
general, the authors estimate that the values of Kg £x identified as 
Experimental Values in Table I, are accurate to within 5% but values of 
K U r are, on the average, probably not better than about 15%. Thus, the 
authors feel that the results for K HEx should be interpreted as trends 
rather than as exact values. In this context, the experimental results 
tabulated in Table I may be sunmarized as follows: 


1 

f a/t 


0.2 




For - 

a/c 

c/r 


1 .0 to 1 .5 
0.5 


K HEx = 

'■ 5 K SEX 


[ 2r/t 

— 

0.5 to 0.9 





' a/t 


0.5 I 




For J 

a/c 

c/r 

: 

1.0 to 1.5 j 
0.5 to 2.0 

f 

K HEx “ 

Sex 


2r/t 

>* 

= 

0.5 to 0.9 J 





r a/t 


0.75 1 




For * 

a/c 

c/r 

— 

2.0 

0.8 to 1 .6 

* 

K HEx = 

°- 5 K SEx 


2r/t 

V 


0.5 to 0.9 j 




Moreover, 

except 

for the shallow flaws, for a/c > 1, K s£x 


These results show that K H£x drops off and K SEx increases as the cracks 
grow deeper. This growth trend is believed to be due to a variation in 

constraint distribution with flaw depth. For the 
parametric ranges included in this study, a/t seems to be the 
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the dominant parameter but a/c is also seen to be important, especially for 
the larger values (i.e., a/c -+• 2.0). For this reason, the tests are 
grouped according to a/t values, and within each group are ordered according 
to increasing values of a/c. 

ANALYTICAL COMPARISONS 

In the middle of Table I, values of the SIF are estimated for each 
case tested using the theories of Hall and Finger, Bowie, and Hsu and 
Liu. (See Appendix A) These results are summarized in abbreviated nor- 
malized form in the Sub-Table at the bottom of Table I for purposes of 
comparison. 

In Group II (a/t - 0.5, a/c - 0.9 to 1.5) all of the theories 
substantially overestimate the SIF both at the surface and at the hole. 
Although the overestimates range from 10% to 220%, they seem to average about 
30% above the experimental result. These same trends are observed in 
Group I (a/t * 0.2, a/c - 1 to 1.5) at the surface and in Group III 
(a/t - 0.75, a/c - 2.0) at the hole. However, for Group I at the hole, 
the Hall-Finger solution underestimates the SIF by about 20%, and for 
Group III at the surface, the Bowie solution under-estimates the SIF. 

As noted earlier, there are many reasons why agreement should not be 
expected between the several theories and experiments. The Bowie solution 
is two dimensional and applies to a very different geometry than the 
cases studied here which are found to be highly geometry dependent. 

The Hall- Finger approach assumes that fracture initiates at the hole 
boundary under a state of plane strain constraint and obtains empirical 
functions from experimental data only for cases where a/c < 1 .0. The 
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Hsu-Liu Theory provides a highly empirical modification of the Bowie 
solution for obtaining the SIF both at the surface and at the hole and 
neglects the back surface of the plate altogether. As noted earlier, 
the effect of increasing a/t was to decrease and to increase K^. 

The decrease of with increase in a/t was also observed by Hall and 
Finger as an "unexpected result." The authors conjecture that this result 
may be due to the fact that, when the crack is deep, a substantial 
part of the load is transferred to the side of the hole opposite the crack 
due to the greater stiffness of this part of the plate. The remaining 
ligament between the crack and the back surface of the plate thus carries 
a reduced load and this effect is in evidence in the results shown here. 

Values of a/c are low enough in tests 1, 3 and 4 that some correlation 
might be expected with the Hall -Finger Criterion. However, differences 
between that criterion and experimental results run as high as ±20%, 
indicating that substantial three dimensional effects may be averaged out 
in the Hall-Finger approach. 

PROBLEM CHARACTERIZATION AND CONCLUSIONS 

It is clear from the results of the Sub-Table under Table I that none 
of the theories studied here are (in their published form) suitable for 
dealing with deep flaws where a/c =>1.0. It is also clear that this problem 
is highly three dimensional and strongly geometry dependent. In order to 
assist the designer in converting these results into an interim design philos- 
phy until further studies can be carried out, the following suggestions are made 
In Reference [1], the crack growth occuring in this type of problem 
through the thickness is estimated from the SIF for a quarter circular 
flaw emanating from the corner of a quarter infinite plate. It is 
clear from the value of « HEX /a c 1/2 from the tests in Group I that such 
a solution should be corrected to account for the stress raiser effect 


9 



due to the hole and that this correction should be of the order of 2.0 
to 2.5. Although difficult to justify on purely analytical grounds, 
the correction obtained from the Bowie solution as c->0 by Hsu and Liu 
appears to be somewhat conservative but not unsuitable. here. For . Groups 
II and III, failure to correct for the presence of the hole is apparently 
more than offset by the load transfer mechanism mentioned earlier. A 
20% or so decrease in the value of C in the Hall-Finger criterion would 
yield fairly accurate Kj estimates for K^. 

Once the crack has grown through the plate thickness, it forms a 
Bowie Crack for which Kj (or K g ) is accurately predicted by the Bowie 
Theory. In fact, the use of the Bowie Theory to predict before the 
crack breaks through to the far side of the plate could be used {if Kg is 
decreased by about 20%) for Groups I and II. For Group III, the Hsu-Liu 
Correction would again be necessary. 

The problem studied here is a complex, three dimensional problem and 
the present study can only be expected to identify trends and point the 
way for future research. The trends and conclusions from the present 
study may be summarized as follows: 

1) Extension of a crack emanating from a corner of intersection of 
a hole with a plate under monotonically increasing load is not self- 
similar. (In this work, most of the growth was through the thickness*) 

As the flaw depth increases decreases and increases. 

9) Existing theories and design criteria significantly overestimate 
the SIF at both the hole and the surface except for shallow flaws at the 
hole and deep flaws at the surface. 

3) As an interim design criterion, the Bowie Theory with magnifications 

of the type proposed by Hsu and Liu for a/t = 0.75 or greater is tentatively 

*Test results on fatigue crack growth at U. S. Air Force Materials Laboratory 
show that, in some cases, significant growth may occur in the plate surface. 
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recommended. For estimates, a Hall-Finger approach with adjusted 
coefficients and relationships is suggested. 

The authors are planning to extend the current study to other crack 
geometries with and without hole fasteners in the near future. 
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APPENDIX A - Related Theoretical Solutions 


Two approximate theories have been proposed for the problem 
described in the foregoing and the two dimensional theory of Bowie 
was also used. These theories will be briefly described here and their 
limitations noted. 

I. The Bowie Solution [3] — Bowie utilized conformal mapping 

in order to obtain a solution to the two dimensional problem of a through 
crack emanating from a hole under remote extension. The mapping function 
was constructed by considering separately the transformation between two 
upper half planes which carries the real axis of one into the real axis 
interrupted by branch cuts of finite length of the other, and the transfor- 
mation which maps each of these half planes into circles and their 
exteriors. The mapping function used was the product transformation and it 
was approximated by a polynomial. The resulting SIF for a single crack 
normal to a field of remote uniaxial tension may be expressed as 

K, = 5/5T F(t) (A-l ) 

where 

a is the remote tension normal to the crack 

L is the crack length 

r is the hole radius 
F(^) is numerically evaluated. 

II. Hall and Finger [2] - - These investigators prepared a series 
of test specimens in accordance with Table A-l. Notches were inserted 
using an EDM and specimens were fatigued to start the cracks. Residual 
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static strength tests were then conducted. Hall and Finger then postulated 
that catastrophic fracture would originate at the intersection of the flaw 
border with the hole and that this region would be under a plane strain 
type of constraint. On this basis they proposed, as a fracture criterion: 

ava> (a,c,r,t) = K Jc (A-2) 

Then, on the basis of parametric crossplots they refined Equation (A-2) 
into the form: 

Co /a F(| r ) • G(f) • H(£) • K, c (A-3) 

where C=l.l, H(jr) « 2/r/t, and F and G were obtained graphically from 
the test data. They obtained agreement with test data to within 10%. 

Two points however, should be noted: 

(1) Since C, F, G and H were obtained from the test data, 

Eq. A-3 is necessarily limited to the range of test data which includes 
no a/c >1.0. 

(2) While the assumption that plane strain constraint exists at 
the point where the flaw intersects the boundary is certainly the safe, 
logical assumption to make for a design criterion, and seems to yield 
reasonable results in their work, the complex variation in Kj and the 
constraint along the flaw boundary is masked by this assumption so that 
one would not be surprised to see deviations from this criterion for 

other geometries. Nevertheless, this study should be useful in establishing 
trends due to variations in test geometry. 

III. Hsu and Liu [5] — These investigators modified the elliptic 
flaw solution to the form: 

K t = a/rTb - * B • M, (A-4) 

| — -r ■ 
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where: 


b = c in the notation of this paper. 

B = factor which corrects crack solution for the presence of the 
hole. It consists of F(L/r) in I. 

Mi = ) sin g + cos g} where M-j is a front surface correction 

factor and g is measured from the front surface to a point on the flaw 
border. 


$ = 


■ f‘i 


r c 2 

( - -y- 5 - ) sin^g j dg 


Unfortunately, no attempt was made to account for the influence of 
the back surface of the plate upon the stress intensity factors. Moreover, 
the analytical basis for the B factor is two dimensional and it's 
variation along the flaw border is quite arbitrary and empirical. 
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TABLE A-l* 


TEST PROGRAM FOR FLAWS ORIGINATING AT HOLES 


■ 

■ 

a 

t 

Number of Te 

sts For a 

/c « 

■ 

0.25 


1.00 



0.2 

2 

2 

2 


1.0 







0.5 

2 

2 

2 



0.8 

2 

2 

2 

-320 








0.2 

2 

2 

2 

2219-T87 






AL. 

0.5 







0.5 

2 

2 

2 



0.8 

2 

2 

2 



0.2 

2 

2 

2 


1.0 

0.5 

2 

2 

2 

-320 


0.8 

2 

2 

2 








0.2 

2 

2 

2 

Ti. -5A1. 







0.5 

0.5 

2 

2 

2 

-2.5Sn 








0.8 

2 

2 

2 


*From Reference [2] 
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Figure 1 - Problem Geometry and Stress Notation 
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2 






(b) 5th multiple 


Figure 4. Typical Fringe Patterns 
(Dark Field, 15X) 
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NORMALIZED APPARENT STRESS INTENSITY FACTOR 
K &p-s / <5‘(c)' /2 OR K ap _ h /<t(ci) i/z 




TABLE I - DATA AND RESULTS 


Group I Group II Group III 


Test Number 

1 

Bl 


4 


H 


Bl 

Ihbhb 

Hole Radius (r) in. 

0.25 

0.25 

0.19 

0.25 


B 

0.25 

0.19 

HR1 

Plate Depth (t) in. 

0.58 

0.93 



0.78 

nm 

KB 

0.78 

m 

Crack Length (c) in. 

0.11 

0.14 

0.38 

0.25 

0.30 

B 

0.20 

0.29 


Crack Depth (a) in. 

0.10 

0.22 

0.33 



0.25 

0.39 

0.60 

B 

Remote Stress (a) psi 

14.38 

15.62 

12.34 



12.78 



7.06 

a/t 

0.18 

0.24 

0.43 

0.48 

0.49 

B 



0.55 

a/c 

0.95 

1.59 

0.86 

1.10 


— m 

Hi 


B 

0.42 

c/r 

0.43 

0.56 

2.02 

0.98 


B 

0.79 

1.56 


2r/t 

1 0 . 87 

0.54 

0.49 

0.88 

0.48 

0.94 

0.95 

0.48 

0.48 

Stress Intensity Factors: 

Experimental Results 



Surface K SEx psi-in 1 ^ 2 

9.43 

12.84 

8.89 



10.00 

12.23 

8.79 

6.89 

Hole K HEx psi-in^ 

13.88 

18.26 

9.31 



Bl 

IB! 

Bl 

IB 

Empirical Results 

Hall-Finger K^p psi-in 1 ^ 2 

10.85 

13.57 

11.30 

10.59 

10.37 

11.92 

IB 


8.23 

Bowie Kg psi -in ^ 

12.13 

13.80 

m 




9.56 



Hsu-Liu (Surf.) psi-in 1 ^ 

12.19 

16.25 

11.53 

11.03 



IB 

B 


Hsu-Liu (Hole) K HLH psi-in 1 ^ 

17.06 




27.80 

19.38 

22.94 

17.64 

15.09 


24.00 

SUB-TABLE 

Test Number 

1 


B 

4 

5 

6 

7 

— 

8 

IEH 

KBPBBBI 

w 5c ’ /2 

2.0 

2.2 

1 .2 

1.7 

1.2 


B 

1.8 



2.6 

2.4 

1.5 

2.0 

1.6 





K HLS /aC ' 

2.6 

2.8 

1.5 


B 





W 5a ' /2 

3.0 

2.5 

1.3 

1.6 

1.0 



KB 

B 

- T/2 
K HF /aa ,/£ 

2.4 

1.9 

1.6 

1.9 

1.4 

1.9 

1.6 

1.2 

1.8 

_ I/? 

W OT 

3.7 

2.7 

4.0 

3.5 

3.2 

2.8 

2.4 

2.4 

5.2 


♦dropped due to data scatter 
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